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Abstract 

The mass splittings of the pseudoscalar and vector D and B hght-heavy 
quark systems have been calculated using the method of QCD sum rules. Elec- 
tromagnetic, quark mass, and nonperturbative QCD effects are all included. 
The results are in good agreement with experiment. A measurement of isospin 
splitting for the vector B mesons would give valuable information about quark 
mass splittings. 
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The origin of mass differences in isospin multiplets has long been of great interest 
in nuclear and particle physics as a source of information about symmetry violations. 
Hadronic isospin violations are particularly important in that they arise from non- 
perturbative Quantum Chromodynamics (QCD) as well as quark mass differences 
(see Ref.[0] for a review of the early work in this area), and of course electromag- 
netic effects. Among the first applications of the method of QCD sum rules was 
the study of isospin violations in the p — uj system 0, in which it was recognized 
that the isospin splitting of the light-quark condensates can produce effects as large 
as the current-quark mass splittings and electromagnetic effects. Recently the QCD 
sum rule method has been used to study the neutron-proton mass difference [|, |[, 
the octet baryon mass splittings and the mass differences in the charmed meson 
systems (the D and D* scalar and vector mesons) 0. 

In the present work we study the mass splittings of the D,D*, B and B* systems 
using the QCD sum rule method. A main feature of the study is a consistent treat- 
ment of electromagnetic corrections within the framework of the method. Rather 
than make use of estimates of electromagnetic effects from quark models as used in 
previous work, or a phenomenological parametrization of electromagnetic effects as in 
Ref , we explicitly calculate the two-loop electromagnetic corrections as well as the 
perturbative and nonperturbative QCD processes including quark masses up to di- 
mension eight. This is part of a study of light-heavy mesons [0] in which the accuracy 
of heavy quark effective theories is tested. 

The isospin splittings are particularly interesting in that the electromagnetic ef- 
fects and quark mass splittings enter with different relative signs in the charmed vs 
the bottom systems, magnifying the role of the isospin violation of the quark conden- 
sates. On the other hand, nonperturbative QCD mechanisms are expected to be less 
important for heavy quark systems, so that an accurate treatment of electromagnetic 
processes allows one to study light quark mass differences. That will be our main 
conclusion. 

One can write the correlators in the QCD sum rule as 

cx 

Yitip') = Uo{p^) + rrigUUp') + -^Uemip^). (1) 
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For no(p^), the leading term for the light heavy quark system, the isospin violation 
comes from the isospin splitting of the quark condensates in the non-perturbative 
power correction, which has been given in previous publications]^, Ij. The second 
term in Eq. |I| is the light quark mass, rriq, expansion of the total correlator; and 
the third term corresponds to the electromagnetic corrections, which have been not 
treated consistently in the framework of the QCD sum rule method in the literature. 
The rriq expansion is extended to dimension 5 in our calculation, and furthermore, 
we also include the one-loop corrections to the quark condensate, which we believe is 
significant. Thus, Umip"^) in Eq. |l] can be written as 

n^(g2) = Cjl + Csiqq) + " G)q) (2) 
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where the coefficient Cj represents the perturbative contributions, whose Feynman 



diagrams are shown in Fig. 1. Based on the results in Ref. 
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and l{x) = — /Q^ln(l ~ u)^ is the Spencer function. Particular attention has been 
paid to the rriq expansion of the quark condensate since we find that it plays an 
important in the heavy-light quark systems. In particular, we have considered the 
one-loop corrections to the quark condensate whose Feynman diagrams are shown in 
Fig. 2-a, b, c, and d, the result to the order a, is 
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for the vector currents. The details of the one loop calculations will be given later []^. 
However, we find that the expansions of the one loop corrections to the quark 
condensate are finite after the mass renormalization. The rriq expansion for the quark 
gluon condensate corresponding to Fig. 2-e and f is 
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for the pseudoscalar and vector current, respectively. 

For the electromagnetic corrections, we write Ue.miQ"^) as 



Dil + D3{qq) + 



(10) 



The coefficient Dj in Eq. |T0| is a two-loop perturbative contribution which replace 
the gluon lines in Fig. 1 by the corresponding photon lines. Therefore, the imaginary 
part of the coefficient Dj can be obtained from the two- loop integral in QCD for 
the charge neutral current, where the charges for the heavy quark eg and the light 
quark are equal with opposite sign: 
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(12) 



X 



for the pseudoscalar and the vector currents, where f{x) is given by Eq. ^. The 
calculation of Dj for the charge currents is more complicated. Since the charges for 
the heavy and light quarks are not the same, the Ward identity can not be used here; 
thus, one should do the wavefunction and mass renormalizations at the same time. 
We can write the imaginary part of the correlator for the charged currents as 



ImDi = -CQCqD^ + (cq + eQeq)DQ + (e^ + eQeq)Dq 
where can be obtained from Eqs. |1^ and |ll| 
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Dq q is the integral for the self energy diagrams for the heavy and the light quarks. 
The corresponding imaginary parts are 
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for the pseudoscalar and vector currents, respectively. 

The in Eq. ^ represents the one-photon-loop corrections to the quark con- 
densate, whose Feynman diagrams are shown in Fig. 2. For the charge-neutral 
systems, with the appropriate change in coupling constants, the results are the same 
as the one-gluon-loop corrections to the quark condensate which have been calculated 
exphcitly0. We find 
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for the pseudoscalar and vector respectively. We set yU = Mq in our numerical evalu- 
ations. There are different expressions for the charged -Ds's. We find 
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for the charged pseudoscalar and vector currents of D mesons, and 
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for the charged pseudoscalar and vector currents of the B mesons. There are also 
corresponding one-photon-loop corrections to the quark-gluon condensates and higher 
dimension power corrections. Whether they are important numerically are remained 
to be studied. As we indicate later, the one photon loop corrections to the quark-gluon 
condensate might be important. 

Defining the quantity lJ^ = Mq — for the Borel transformation, where u;^ 
measures the off-shell effects of the heavy quarks, we find 
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for the pseudoscalar and vector current. The correlator Iit{oj%) can be divided into 
the perturbative and nonperturbative parts: 
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where the perturbative term is written in dispersion integral form, and the nonper- 
turbative correlator n„p(cij^) comes from the standard Borel transformations of the 



terms associated with the condensates. The masses Mp and My in Eqs. ^ and |26 
can be obtained from 



M' = Ml^ 00% ^^^H^- (28) 

The mass differences between different isospin states are obtained by adjusting the 
parameters Uq in Eq. to insure that the mass difference is independent of Borel 
parameter uj%. The range of u"^ is chosen in such way so that the the continuum and 
the nonperturbative contributions are minimum. 

There are two parameters that are important for the isospin mass splittings: the 
up and down quark mass difference, rud — mu, and the isospin splitting of the quark 
condensate, which is defined by the parameter 7; 

7 = j^^-l. (29) 

The quark mass difference has been estimated via broken chiral symmetry IH], while the 
parameter 7 in various estimates ranges from — 0.0079[Q to —0.0020, |^. Physically, 
the sensitivity of the isospin mass splittings to the parameter 7 measures the strength 
of the nonperturbative effects on the light-heavy quark system. To examine the 
sensitivity of the mass splittings on the parameters 7, we calculated the mass splittings 
for different 7's in our framework, with the quark mass difference — m„ fixed at 
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3.8 MeV. Our result for different 7 and — m„ fits approximately to the empirical 
formula [all mass differences are given in MeV] 

Mi5± -M^o = 0.1-507 + l.l(md-m„), (30) 

Md*± - Md,o = -0.9 - IO7 + 1.0(md - niu), (31) 
for D and D* mesons and 



Mb± - Mbo = 1.45 - 1.16(md - m„) (32) 

Mb*± - Mb-o = 1.36 - 1.28(md - m„) (33) 

for B and B* mesons. Comparing to the corresponding formula for the neutron-proton 
mass difference 0], neglecting electromagnetic corrections, 

M„ -Mp = -4167 + 0.197(md-m„), (34) 



we find that the isospin mass splittings are much less sensitive to the 7 for the 
D mesons, and the mass splittings for the B systems is almost independent of 7. 
This result differs from the conclusion of Eletsky and loffe who found a stronger 
dependence of isospin mass splittings for the D mesons on the parameter 7. The 
difference is due to their neglect of higher dimension terms, especially the quark 
gluon condensate in no(p^) and Umip^) of Eq. 1, which is very important in our 
calculations. The weak dependence of the isospin splittings on the isospin splitting of 
the quark condensate is consistent with conclusions of heavy quark symmetry that the 
nonperturbative effects become less important for heavier quark systems. Therefore, 
the heavy-light quark systems provide an excellent testing ground for the up and 
down quark mass difference if electromagnetic corrections are calculated accurately. 
For 7 = —0.0079 and rrid — ruu = 3.8 MeV, we have 

M± - M° - I ^-^^ ° (35) 



and 



This is in excellent agreement with data (M/^i — M^o = 4.7 ± 0.7 MeV and M^j.i — 
Md.o = 2.9 ± 1.3 MeV|T^) for D mesons, while the magnitude of our result for 



B mesons is quite large comparing to the data —0.1 ± 0.8 MeV|TO|. One of the 
important features of our calculation is the differential vector-pseudoscalar isospin 
mass splittings A(M); 



A(M) = (M*± - M*°) - (M± - M"), (37) 
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which provides an very important probe of the hyperfine sphttings of the heavy-hght 
quark potential Our result indicates that A(M) is less dependent on the quark mass 
difference rrid — niu, and from Eqs. ^ and we have 

^^M) = 1 ^ mesons 

1 —0.5 for B mesons 

for rrid — rriu = 3.8 MeV. This is in good agreement with recent studies by Cutkosky 
and Geiger|TT]. In particular, the A(M) for B mesons is negative for positive uid—m,,. 



which is opposite to the quark model predictions and several other investigations [p!^ 
A measurement of the isospin mass splittings for the vector B mesons will provide 
a crucial test in this regard. A study in chiral perturbation theory finds |T^ that 
the contributions from quark mass difference only gives 0.3 MeV for D mesons, and 
the rest might comes from the electromagnetic corrections. Our calculation is con- 
sistent with this conclusion, and furthermore we also find that there is only a small 
contribution from the difference of the quark condensates to A(M). 

We find the electromagnetic corrections to the isospin splittings to D and D* 
mesons are quite small. They contribute less than 1 MeV. As a result, our splittings 
in the charmed systems are too small in comparison with experiment for parameters 
giving good fits to the B splittings. One source of this discrepancy might be that 
the one photon loop corrections to the quark gluon condensate are needed, since we 
find that the quark gluon condensate is as important as the quark condensate in 
evaluating the masses of the heavy ligth systems. This might both help to explain 
the electromagnetic corrections are small for the D mesons and improve our fit to the 
B mass splitting with m^ — m„ about 3.0 MeV. This is a somewhat subtle calculation, 
which we do not attempt here. 

In summary, we show how to include the electromagnetic corrections in the QCD 
sum rule analysis for the light-heavy quark systems. This approach can also be 
extended to the study of the neutron-proton mass difference as well as other isospin 
mass splittings. Our results show that the mass splittings are not sensitive to the 
isospin splitting of the quark condensate, and thus heavy-light quark systems are an 
ideal place to study the up and down quark mass differences. We expect the method 
will give a very good qualitative even quantitative descriptions of the isospin-splittings 
for the light-heavy quark systems. Our result for A(M) is in good agreement with 
recent phenomenological studies by Cutkosky and Geiger. Further investigation of the 
electromagnetic effects of higher dimensions terms is needed, and this is in progress. 
Our investigation shows that the isospin mass splittings of the light heavy quark 
systems provide an very important opportunity to study the light quark masses and 
the hyperfine interactions of the quark-quark potential. 

This work is supported in part by National Science Foundation grant PHY- 
9023586. 
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Figure Captions 

1. Processes for two loops correction. 

2. Contribution from nonperturbative processes. 
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